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1. SCOPE 


1.1 GENERAL 

This document is the "as-built" design specification of the 
CAMS/CAS Interface Tape Report Generation Program for LAC IE 7. 




2. APPLICABLE DOCUMENTS 


• TIRF 77-0052 

• TIRF 77-0040 

• Specification for the CAMS/CAS Interface Tape Report 
Generation Program - LEC-9151 

• CAMS/CAS Interface Control Tape format specification in 
Earth Resources Data Format Control Book Volume 1 (PHO-TR543, 
Rec. A. Change 3) 

• TIRF 78-0010 

• CAMS/CAS Interface Tape Printout after LACIE 6A JSC Memo- 
randum, SF4-77-7-13 , 7/21/77. 

• "As-Built" Design Specification of the CAMS/CAS interface to 
tape Report Generation Program LEC-11292 

• Acceptance Test Specification For CAMS/CAS Interface tape 
report generation Program for LACIE 7-LEC 11787 

• Clarification and Prioritization of LACIE 7 CCIT Report JSC 
memorandum, SF4-77-11-8, 11/4/77. 


3. SYSTEM DESCRIPTION 


3.1 HARDWARE DESCRIPTION 
N/A 


3.2 SOFTWARE DESCRIPTION 

The purpose of this program is to produce CAMS reports from 
data on the CAMS/CAS interface tape. 

3.2.1 SOFTWARE COMPONENT NO. 1 (CAMRPT) 

The main program CAMRPT reads control cards, locates segment data 
on the input tape and calls subroutines to generate requested 
reports. 


3. 2. 1.1 Linkages 

CAMRPT calls subroutines CDRED, CAMREC, BAUEXT , CLURES , STDATA , 
CONRED, and DOTRAY. 

3. 2. 1.2 Interfaces 
N/A 


3.2. 1.3 Inputs 

CAMRPT control cards are: SEGMENT XXXX, RECORD ID XXXXXX 

XXXXXX, ALL, END. CAMS/CAS interface tape records are inputs to 
CAMRPT. See reference 4 in section 2, for record formats. 

3.2. 1.4 Outputs 

An error message is output indicating a bad data card. If a 
requested segment is not on the input tape, the program writes 
a message to that effect. 
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3. 2. 1.5 Storage Requirements 
Total space allocated is 2390 bytes. 

3.2. 1.6 Description 

CAMRPT is the CAMRPT main program. The program sets the printout 
option indicator PRTOUT to 0 initially. In this mode the output 
of some reports is conditional. On the first call to tape read 
subroutine CORED, the program reads two data cards specifying the 
device code (M or X) and the unit number (0 to 1), then calls 
CONRED to read constants to be used in the calculations for 
DW & DS. CAMPRT next reads a program control card and tests the 
first non-blank character for one of the following: S, R, A, or E. 

If the card is blank or the first character is not one of the 
above, the program prints an error message on the line printer 
•and stops. 

The action taken for each control card is given below. Note 
that if a control card other than E is read in, the printout 
option indicator PRTOUT is set to 1. This is the option to out- 
put all reports, including conditional reports. 

e S - Option indicator PRTOUT is set to 1. The program obtains 

the segment number from the input card. The program searches 
the input tape for a recognition segment record whose 
segment number matches the control card segment number. 

If a match is not found, a message is printed and the 
program goes to read the next control card. If a match is 
found, the reports for the segment requested are 
generated. The program then reads the next control card. 

• R - The action taken is the same as for the S card above, 
except that the record identification number is used 
instead of the segment number. 

e A - Option indicator PRTOUT is set to 1. Beginning with the 
segment on the tape at which the tape is currently 
positioned, the progran generates reports for that 








segment and all the following segments. When the second 
tape end of file, indicating end of data, is roached, the 
program rewinds the tape and returns to read the next 
control card. 

e E - If the printout option indicator PRTOUT is 0, the program 
generates reports in the limited printout mode, rewinds 
the tape, and then stops. If PRTOUT is 1, the tape is 
rewound and the program stops. 

The tape read subroutine CDRED tests all records for valid 
characters. If any invalid data is Encountered an error messaae 
is output to the line printer to inform the user that the data 
for that segment or record is questionable. Any invalid characters 
in the record are converted to ones (1) and normal processing 
is resumed. 

To generate reports for a segment the program first calls CAMREC 
to process recognition segment records. The processing entails 
generation of the optional classification Data report, the standard 
Separability Report and the standard header sheet for the report. 

After CAMREC, the programming calls BAUEXT to process the Bias 
correction results records. Data from the Bias correction results 
records is saved in coinpton blocks BIAS and Dummy. CAMRPT next 
calls CLURES to process the clustering results records and to 
generate the conditional cluster report. In addition CLURES 
saves cluster dot data in common block CLDOT for later use in the 
optional Dot report. Next STDATA is called to process the statistics 
records and generate the optional Statistics Report. Finally 
DOTRAY is called to process the Dot subset records and generate 
the standard Dot Label/Classification, Bias Correction Classifies** 
tion, Dot Label/Cluster and Bias Correction Cluster Reports. 

In the limited printout mode the decision as to whether or not 
to output the optional reports is not made until the. percentage of 



correctly classified Dots (PCC-1 & PCC-2 for Bias Correction 
Classification report) or percentage of correctly clustered Dots 
(PCC-1 & PCC-2 for Bias Correction Cluster Report) are calculated 
in subroutine DOTRAY. If any of the values are less than 80%, 
PRTOUT is set to 2 in DOTRAY and the conditional reports are read 
from disc and written to the line printer. In the full printout 
mode (PRTOUT-1) the conditional reports are always retrieved from 
disc and printed out. If PRTOUT was ■ 2 it is reset to 0 after the 
report for a segment has been butput. 

3.2. 1.7 Flowcharts 
See Flow Diagram 1. 


3. 2. 1.8 Listing 
See Appendix A. 


3.2.2 SOFTWARE COMPONENT NO. 2 (CAMREC) 


This program procaaaass classification results contained in 
recognition segment records, outputs the conditional Classifi- 
cation Data report, saves data from the recognition segment . 
records for use in generating the standard Separability Report 
and the standard CAMS Interface Report Header sheet. 

3. 2. 2.1 Linkages 

CAMREC is called by CAMPRT and calls subroutines BIAPRT, CPIPO, 
MV, CORED, BN? and SEPRPT. 

3. 2. 2. 2 Interfaces 
N/A 

4 

3. 2. 2. 3 Inputs 

Recognition segment records, containing subclass a priori and 
threshold values, subclass related classification results, and 
feature selection Bhattacharyya separability data for available 
acquisitions. 

3 . 2 . 2 . 4 Outpu ts 

CAMS Interf o .e Report Header and a Classification Data report. 


3.2. 2.5 Storage Requirements 
Total space allocated is 2845 bytes. 

3. 2. 2. 6 Description 

CAMREC is called with the first recognition segment record for 
the segment to be processed residing in array IBUF. CAMREC 
first calls BIAPRT with PASS-1 to have the report heading, segment 



number, record ID, and acquisition dates output. Next CAMREC 
saves the number of channels used in classification and the 
Bhattacharyya separability data from the first recognition seg- 
ment record for later use in generation of the Separability Report. 

Title and column headings for the classification section of the 
report are written out by CAMREC. Processing of classification 
results begins by setting the location in array IBUF of the 
first subfield containing subclass related results. Subfield 
contents are accessed by calling CPIPO. CPIPO returns the class 
portion of the subclass name and the counts PI and PO of pixels 
classified into, and thresholded out of the subclass. If the 
first character of the class name is X, PI is added to the X 
category pixel count. If the category is W, for wheat, then the 
count for the first wheat class is set to PI and the wheat class 
name is saved in CLIST. PO is added to the total of pixels 
threshold, TC, in the COMMON blocks CBIAS. 

In processing for the second, and subsequent subclasses, the 
program calls CPIPO to get the next class name, checks to see 
if it is wheat, and, if so, compares it to the last class name 
in CLIST. If it is not the same, the new name is saved in 
CLIST and the class index is incremented by 1. This causes 
wheat class pixel count PI to be tallied in the next results 
array location. 

After all classification data has been processed the feature 
selection Bhattacharyya separability data is saved from the last 
recognition segment record and SEPRPT is called to generate 
the normal Separability Report. 

3. 2. 2. 7 Flowcharts 
N/A 


3. 2. 2. 8 Listing 



See Appendix A. 
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3.2.3 SOFTWARE COMPONENT NO. 3 (BAUEXT) 

This program saves data contained in the clustering bias 
correction and classification bias correction results records. 

3. 2. 3.1 Linkages 

BAUEXT is called by CAMRPT. It calls CORED. 

3.2. 3.2 Interfaces 
N/A 

3.2. 3.3 Inputs 

Clustering Bias correction and Classification Bias Correction 
results records. 

3.2. 3.4 Outputs 
None. 

3.2. 3.5 Storage Requirements 
Total space allocated is 2306 bytes. 

3.2. 3.6 Description 

BAUEXT is called when the main program reads the first clustering 
Bias correction result record. BAUEXT saves the following data 
from both the Clustering Bias Correct in and Classification Bias 
Correction results records for up to 8 categories of interest 
plus the "designated other” and "unclassified" category" 

Pixel pouplation 
Bias corrected estimator 
Machine estimate 
Random estimate 

Variable of bias corrected estimate 
Variance 



In addition the number of categories of interest and the character 
used for the categories of interest are saved. All data is saved 
in common blocks CBIAS arrays. This data is used by CALC for 
certain calculations and by BIAPRT for output of the normal 
Bias Correction reports. 

3.2. 3.7 Flowcharts 
N/A 

3.2. 3. 8 Listing 
See Appendix A. 











3.2.4 SOFTWARE COMPONENT NO. 4 (CLURES) 

This program process the cluster results records and generates 
the conditional cluster report. 

3. 2. 4.1 Linkages 

CLURES is called by CAMRPT. It calls CORED. 

3.2. 4.2 Interfaces 
N/A 


3.2. 4. 3 Inputs 
Cluster results records. 

3. 2. 4. 4 Outputs 

The conditional Cluster report. 

3.2.4. 5 Storage Requirements 
Total space allocated is 5632 bytes. 

3. 2. 4. 6 Description 

CLURES is called when the main routine reads the first cluster 
results record. The program decodes ALSETS, the total number 
of clusters, and SETSR, the number of clusters in the current 
record. The routine then outputs the cluster report header, 

ALSETS as clusters generated and any options used. Next CLURES 
saves all data for each cluster for later output. When all 
clusters in the current record have been processed, another cluster 
results record is read in and processed as above. 

After all cluster have been processed and if the cluster/dot 
report option is set, four additional cluster results records 
need to be processed. The processing .consists of saving all dot 



information in an array called DOTBUF for later output. In 

addition each dots cluster assignment is transfered to the 

common blocks CLCOM for later use in a different report. Finally 

when all cluster results records are processed in the above 

manner the cluster information is output as follows. For each 

cluster the program outputs the cluster name, the Labeling dot 

match name, Ll distance, categorie used , brightness ~afid 

* ... 

greenness numbers for all Acquisitions used and information on 

all dots in the cluster. The clustering channel list is 

written at the end of the report. 

3. 2. 4. 7 Flowchart 
See Flow Diagram 2. 

3. 2. 4. 8 Listing 
See Appendix A. 



3.2.5 SOFTWARE COMPONENT NO. 5 (STDATA) 

This subroutine formats and. outputs field and subclass statistics 
data. 

3.2. 5.1 Linkages 

STDATA is called by CAMRPT. STDATA calls subroutines KNT, MDTTL , 
MEAN, POP, CDRED, SNAME, FANME, STDMP, and BNT. 

3.2. 5.2 Interfaces 
N/A 


3. 2. 5. 3 Inputs 

The statistics record, containing, for fields or for subclasses, 
the population and values of the mean and standard deviation 
by channel* 

3. 2. 5. 4 Outputs 

The conditional statistics report. 

3. 2. 5. 5 Storage Requirements 
Space allocated is 3792 bytes. 

3. 2. 5. 6 Description 

STDATA is called from CAMRPT. By means of decode statements, 
the program converts several variables from input character 
format in IBUF to integers. The variables are ALSETS, the total 
number of statistics sets, SETSR, the number of sets in the 
current record, and NCH, the number of channels. STDATA calls 
souboutines to move data from input record subfields to print 
buffers. SNAME and FNAME move name data and insert SUBCL and 
FIELD designations in the print buffer. POP is called to move 









population data. NDTTL is called to supply column headings for 
means and standard deviations, which are transferred to a print 
buffer by MEAN. MEAN also puts decimal points where needed. 

The variable DSETS, set to 5, controls the number of statistics 
sets to be accumulated before outputting the print buffers. 

When the current record statistics sets counter reaches SETSR, 
and ALSETS sets have not yet been processed, STDATA calls CDRED 
to read the next statiscics record from tape. In addition, the 
pixel population and classified percentage for the category 
"unassigned" in the bias correction cluster report is calculated 
and saved in the common blocks CBIAS. 

3. 2. 5. 7 Flowcharts 
N/A 


3. 2. 5. 8 Listing 
See Appendix A. 



3.2.6 SOFTWARE COMPONENT NUMBER 6 (DOTRAY) 

This program processes Dot Data records and generator the 
conditional Dot Report. 

3. 2. 6.1 Linkage 

DOTRAY is called by CAMRPT. It calls BIAPRT, CALC and CDRED. 

3. 2. 6. 2 Interfaces 
Dot Data records. 

3.2.6. 3 Inputs 
N/A 

3. 2. 6. 4 Outputs 

The Conditional Dot report. 

3. 2. 6. 5 Storage Requirements 
Total space allocated is 2306 Bytes. 

3. 2. 6. 6 Description 

DOTRAY is called by CAMPRT after reading the first Dot Data 
record. DOTRAY next outputs the report header to the top of 
the next page and processess the Dot Data records until all 
209 dots are processed. DOTRAY also saves the dot Lab le, dot 
type and classification for each dot in the common block CBIAS 
for later use by BIAPRT & CALC. 

The processing involves outputting the following for each of 
209 dots: 

Dot number, line and pixel number for the dot, type and label 
(if any) for the dot, cluster and classification is well as the 
greenness and brightness of up to 4 acquisition for the Dot. 
Each record contains data for 15 dots. After the 15 dots in the 
record has been processed and output the next record is read in 
via CDRED and processed. After all dots are processed DOTRAY 




calls CALC to perform calculation for the bias correction reports 
then calls BIAPRT to output the bias correction reports. 

Finally, DOTRAY check the value PRTOUT. If PRTOUT is 1 DOTRAY 
returns to the main program. If PRTOUT is'not 1 all valid PCC 
values (percentages of TYPE 1 & TYPE 2 correctly classified or 
clustered dots) are tested. If any of the valid PCC's are 
less them 80%, then PRTOUT is set to 2, to indicate to the main 
progreun that the conditional reports are to be read from the 
disk and output to the line printer. 

3. 2. 6. 7 Flowcharts 
N/A 

3.2.6. 8 Listing 
See Appendix A 





3.2.7 SOFTWARE COMPONENT 7 (SEPRPT) 

This program processess the separability data which was saved 
for it by CAMREC and outputs the separability report. 

3.2. 7.1 Linkage 

SEPRPT is called by CAMREC. 

3.2. 7.2 Interfaces 
N/A 

3.2. 7. 3 Inputs 

Segment Recognition record number 1. 

3. 2. 7.4 Outputs 

The normal Separability report. 

3. 2. 7. 5 Storage Requirements 
Total space allocated in 588 bytes. 

3.2. 7.6 Description 

SEPRPT is called by CAMREC after all separability data has been 
saved. SEPRPT then outputs the separability report which 
consists of channel combinations for up 4 
(16 channels) and selection Battacharyya separability data. 

3. 2. 7. 7 Flowcharts 
N/A 

3.2. 7.8 Listings 


See Appendix A 
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3.2.8 SOFTWARE COMPONENT NO. 8 (CALC) 

This subroutine calculates data necessary for the output of the 
Bias Correction Classification and Bias Correction Cluster 
reports . 

3.2.8. 1 Linkage* 

CALC is called by DOTRAY. 


3. 2. 8. 2 Interface 
N/A 

3. 2. 8. 3 Inputs 

Dot Data from common blocks CBIAS. 


3.2. 8.4 Output 
None 

3. 2. 8. 5 Storage 

Total space allocated is 4031 bytes. 


3. 2. 8. 6 Description 

CALC performs 2 identical calculations on different sets of data. 
The first data set involves Classification data and the second 
involves Cluster data. CALC first sets all needed variables 
to zero. If the flag CLADUM is equal to 1 no classification 
calculations are performed and CALC goes directly to the 
cluster calculations. Otherwise CALC checks "category 8 interest" 
for the presence of an N S", "N", or "W", and sets pointers and 
indication as necessary. Next certain arrays are zeroed. 


Following is the variables calculated for each of 209 dots. It is 
understood in every case that the conditions apply to DOTS which 




are labeled, that is, LBLED (I) la not blank, and to DOTS not 
classified as DU or DO. 

• NTYP1 - The number of DOTS which are either type 1 or type 3. 

• NTYP 2 - The number of type 2 DOTS 

• NAII - The number of type 1 type 3 DOTS whose label and 

classification are the same. 

• NGIJ - The number of type 1 and type 3 "Grain TYPE" DOTS 

whose label and classification are not the same. 

• NOCL - The number of type 2 DOTS which are both labeled 

and classified. 

In addition when calculating the above for all dots a bias correction 
vector table is calculated. This table consists of a two 
dimensional array and contains summations of dots that have 
valid labels as the first index and valid classification as 
the second index. 

In addition the following variables are calculated: It is under- 

stood that classified does not mean "threshold". 

• ALGT - The number of TYPE 2 dots labeled either "W" or "S" 

and classified either "W" or "S". 

• ALBG - The number of TYPE 2 dots labeled with any "category" 

used" and classified either "W" or "S". 

e ALNT - The number of TYPE 2 dots labeled in any "category used" 
other than "W" or "S" and classified in any "category 
used" other than "W" or "S". 

• ALNT - The number TYPE 2 dots labeled in any "category used" 

and classified in any "category used" other than 
"W" or "S". 

The subroutines computes the corrected percentages, uncorrected 

populations, variances, uncorrected percentages and random 

/ 

yrf* 'f 



■ample data for all "categories used" plus "grain category". 
The equations for these calculations are shown in Attachment B. 


CALC now performs the identical calculations on the cluster data 
i\£i.nT checking CLUDUM as outlined above. 

3. 2. 8. 7 Flowchart 

See Flow diagram No. 3. 

3.2.8. 8 Listing 
See Appendix A. 



3.2.9 SOFTWARE COMPONENT NO. 9 (BXAPRT) 

This program outputs the report header sheet# Label/Classification 
table# Lable/Cluster table and the Bias Correction reports. 

3.2.9. 1 Linkage 

BXAPRT is called by CAMREC 4 DOTRAY. It calls BNT. 

3. 2. 9. 2 Interface 
N/A 

3.2.9. 3 Inputs 

Segment recognition records and common blocks CBXAS# Dummy and 
CLCOM. 


3.2. 9.4 Outputs 

Report header sheet^TYPEl and TYPE2 Dot Label/Classification 
report# Bias Correction Classification Report, TYPE1 and TYFE2 
Dot Label/Cluster report and the Bias Correction Cluster reports. 

3. 2. 9. 5 Storage Requirements 
Total space allocated in 3712 bytes. 

3.2. 9.6 Description 

When CAMREC calls BXAPRT it sets PASS*1. This causes BXAPRT 
to output the report header which contains tape number# DPAR 
No.# record XD# segment number and all acquisition dates. This 
information is retrieved from the segment recognition record 
which had been read into XBUF. 

When DOTRAY calls BXAPRT it sets PASS-2 . In this mode of opera- 
tion up to 2 similar reports can be generated. Xf CLADUM is 
equal to 1 no classification report is generated. 




Otherwise DOTRAY generates an 11 by 19 matrix of user label/clas- 
sification entries for type 1 and type 3 DOTS. A similar Matrix 
is generated for type 2 DOTS which also includes type 0 DOTS. 

DOTS with a classification label of DU or DO do not appear in 
either matrix. 

Next the bias correction classification report is output using 
data calculated by CALC and stored in common blocks CBIAS. The 
report consists entries for all "categories of interest", DO, 

TH, UN and grain as follows: 

Pixel population, classified and corrected percentages, variance 
and random sample estimate:- , 

Also the alpha value matrix is output, and PCC values, DW and DS 
and lastly the Bias correction vectors and totals. 

A similar report is generated for the cluster data if CLUDUM is 
not equal to 1. 

3. 2. 9. 7 Flowcharts 
See flow diagram 4. 

3.2.9. 8 Listing 
See Appendix A. 
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THE LOGIC FOR SECOND PASS IS THE SAME AS THE FLOW CHART 
ONLY THE VARIABLES ARE CHANGED. 
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4. OPERATING PROCEDURE 


4.1 GENERAL 

CAMRPT is an RSX-11D Batch program which reads a CAMS/CAS Inter- 
face tape and generates a series of reports which are output on 
the line printer. It requires data card inputs. 

4.2 DECK SET UP 

The first 22 data cards define the input tape drive , input tape 
unit and curve fitting constants for DW 6 DS calculations. The 
DS & DS formular are shown in Appendix B. The card formats 
are: 

M or X (tape drive) 

0 or 1 (tape unit) 

10 DW curve fitting constant cards as described in Appendix Cl 
10 DS curve fitting constant cards as described in Appendix C2 

Entries always start in column 1. To execute the CAMRPT default 
option for a limited printout of reports, an END card must 
follow the two data cards above. If the option for a full 
output of all reports is desired, the control card sequence is: 

A (for all reports) 

END 

If the user desires to obtain the output for only a single seg- 
ment on the input tape, segment 9681 for example, the control 
card sequence is: 

S 9681 

END 

To obtain the output for segment 9681, and all segments following 
9681, the control card sequence is: 

S 9681 


The above option is used when there is a bad segment on the input 
tape, to obtain the output for segments following the bad 
segment. 

The Batch deck set up for the CAMRPT default option using input 
tape unit MTO is as follows: 

$JOB/NAME-CAMRPT/MCR/LIMIT«99/ACCOUNT«5050 

$DATA 

N 

0 

10 DW constant cards 
10 DS constant cards 
END 
$EOD 

$MCR REM RSXBAT 
$RUN CAMRPT 
$EOJ 

To run the program, mount the CAMS/CAS Interface tape and enter 
a mount message. 

For MTO the message would be: 

MCR > MOU MTO : /CHA* [FOR] (CR) 

Then load the card reader with the CAMRPT Batch deck and enter 
BAT CR: , to read in the deck. 
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The values Dw and Ds are computed as follows* 

rl r* r* 

Dw * 4 p»l q-1 RpqWpWq 4 “pal TpWp or 0 if Category W not used 

r* j* r* 

Ds • C 2 + p-1 ' g“l UpqSpSq 4 p«l VpSp or 0 if Category 8 not used 

• • i i i • 

1 where (Rpq)p-l q-1, {Tp}p-1, {Opq}p-l q-1, and (^p)p-l are 

constants* supplied by user __ 

- Bias corrected estimate for Category W 
8^ - Bias corrected estimate for Category 8 
f ?2 ■ Machine estimate for Category W 

5 2 ■ Machine estimate for Category 8 
Wj - Random estimate for Category W 

5 3 - Random estimate for Category S 

Vj - Variance of Bias corrected estimate of W 

Sj « Variance of Bias corrected estimate of 8 
WtW 

w 5 “ «:W 4 fi,W 4 -HTU 

e 8:8 

®5 W:8 4 8:S 4 Ni4 

Mfi - S 6 - PCC1 

«> . NiM 

*7 / WiN 4 S’iN 4 tfrg 

M 8 -8 8 - PCC2 
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a. Constants* for Do calculations 

Card 1 contains as follows : 

blanks or XXX.XXXXX in first 10 columns only 

Card 2 contains thru Tg as follows: 

blanks or XXX.XXXXX for aach entry. A m ax im u m of 10 colunns 
wide starting in columns 1# 11, 21, ate. 

Card 3 thru 10 contain to Rg g. R^ ^ thru Rj^g on card 3# 
* 2 ( j thru R^g on card 4, ate. Format same as card 2. 

b. Constants* for Ds calculations 

Card 1 contains Cg as follows: 

blanks or XXX.XXXXX in first 10 columns only 

Card 2 contains thru Vg as follows: 

blanks or XXX.XXXXX for aach entry. A maximum of 10 columns 
wide starting in columns 1, 11, 21, etc. 

Card 3 thru 10 contain ^ 1 °8 8* °1 1 t * ura °1,8 0,1 card 3 ' 

0 21 thru U 2 'g on card 4, etc. Format same as card 2. 

*List of test cards on the following page. 
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